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The Cosmic Microwave Background in an Inhomogeneous Universe 

- why void models of dark energy are only weakly constrained by the CMB - 
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The dimming of Type la supernovae could be the result of Hubble-scale inhomogeneity in the 
matter and spatial curvature, rather than signaling the presence of a dark energy component. A key 
challenge for such models is to fit the detailed spectrum of the cosmic microwave background (CMB). 
We present a detailed discussion of the small-scale CMB in an inhomogeneous universe, focusing on 
spherically symmetric 'void' models. We allow for the dynamical effects of radiation while analyzing 
the problem, in contrast to other work which inadvertently fine tunes its spatial profile. This is 
a surprisingly important effect and we reach substantially different conclusions. Models which are 
open at CMB distances fit the CMB power spectrum without fine tuning; these models also fit the 
supernovae and local Hubble rate data which favour a high expansion rate. Asymptotically flat 
models may fit the CMB, but require some extra assumptions. We argue that a full treatment of 
the radiation in these models is necessary if we are to understand the correct constraints from the 
CMB, as well as other observations which rely on it, such as spectral distortions of the black body 
spectrum, the kinematic Sunyaev-Zeldovich effect or the Baryon Acoustic Oscillations. 



I. INTRODUCTION 



The problem of understanding the physical origin and 
value of the cosmological constant is leading us to recon- 
sider some of the foundational aspects of cosmological 
model building more carefully pp. In particular, it is an 
important fact that, at the moment, the spatial homo- 
geneity of the universe on Gpc scales exists by assump- 
tion, and is not yet an observationally proven fact. Given 
this uncertainty, so-called void models can explain the 
observed distance modulus utilizing a spatially varying 
energy density, Hubble rate and curvature on Gpc scales, 
without any unusual physical fields at late times ffi-R)^] 1 . 
Although introduced before the SNIa data [2], they are 
currently attracting attention as, in some respects, one 
of the most conservative explanations for the dark en- 
ergy problem. These models require an as yet unknown 
mechanism for their formation, probably requiring some- 
thing unusual at the start of inflation to create a model 
close to spherically symmetric on Hubble scales. They 
also are in (dire) need of an explanation for the anti- 
Copernican fine tuning that exists: we have to be within 
tens of Mpc of the centre of spherical symmetry of the 



background [IS] [25J H3] > which implies a coincidenc e of , 
roughly, (40 Mpc/15 Gpc) 3 - 1CT 8 (though see Sec.[lV|. 



1 We have in mind the 'very big void models' which vary gently 
over Gpc scales (e.g., 2 13, 48 ). The density profile reaches full 
width at half maximum around the Hubble scale. Other 'void 
models' - or Hubble bubble models - for dark energy suggest 
that because very big voids of 100's Mpc across are observed 
to exist, we could live in the centre of one of these, giving the 
necessary jump in the distance modulus to fit the SNIa — see, 
e.g., gl H| [25] [28] [39]. These can be tested with sufficiently 
many SNIa 1371 , and are often discussed in a perturbed FLRW 
context. In essence, our analysis here applies to these models too, 
but we don't consider these specifically (though see Sec. Ill CI. 



Nevertheless, given that even worse temporal fine-tuning 
exists in our present understanding of the value of the 
cosmological constant and the coincidence problem, these 
models should be taken seriously despite this sinister 
drawback. 

Observationally, it is very hard to disentangle time evo- 
lution from radial variation for a central observer (which 
we consider here). To make matters worse, our current 
lack of a plausible formation mechanism means that ob- 
servationally, at this stage, we should think in terms of 
directly reconstructing the void, rather than constraining 
a model parametrically. In this sense, we have to consider 
the physical conditions for the centre of the void and 
asymptotically as essentially independent to be deter- 
mined observationally (perhaps assuming that all quan- 
tities have the same spatial profile). A common theme is 
to imagine a void in an Einstein-de Sitter model to make 
things simpler (e.g., [40]); but they are so big as to be 
impossible to create from scale-invariant Gaussian initial 
conditions. There seems no real reason to be restricted to 
asymptotically flat models if we know slow-roll inflation 
has to be reworked anyway. 

What observations can reconstruct the void? SNIa and 
other local observations give us a handle on the curvature 
locally which tells us the Hubble rate and total matter 
density locally, as well as their spatial profiles out to a 
redshift of 0(1). CMB observations help constrain the 
void profile out to large distances, because the CMB gives 
us a very precise measurement of the area distance at z ~ 
1100. Moreover, as we discuss in this paper, the CMB 
tells us the baryon fraction and the baryon-to-photon 
ratio at comoving scales ~ 13 Gpc away from us. To 
fully constrain the possible degrees of freedom, we need 
other observations to tell us how the baryon fraction, 



the radiation energy density and the primordial power 
spectrum vary with radius. 

The physics of decoupling and line-of-sight effects con- 
tribute differently to the CMB, and have different depen- 
dency on the cosmological model. In sophisticated inho- 
mogeneous models both pre- and post-decoupling effects 
will play a role, but Hubble-scale void models allow an 
important simplification for calculating the moderate to 
high I part of the CMB gQl [42] (see also [68]). The co- 
moving scale of the voids which closely mimic the ACDM 
distance modulus are typically O(Gpc). The physical 
size of the sound horizon, which sets the largest scale 
seen in the pre-decoupling part of the power spectrum, 
is around 150 Mpc in comoving units. This implies that 
in any causally connected patch of the Universe prior to 
decoupling, the density gradient is very small. Further- 
more, the comoving radius of decoupling is larger than 
10 Gpc, on which scale the gradient of the void profile 
is small anyway (by assumption of course - it could be 
set to zero). For example, at decoupling the total frac- 
tional difference in energy density between the centre of 
the void and the asymptotic region is around 10% [55] : 
hence, across a causal patch we expect a maximum I % 
change in the energy density in the radial direction, and 
much less at the radius of the CMB that we observe for a 
Gaussian profile. This suggests that before decoupling on 
small scales we can model the universe in disconnected 
FLRW shells at different radii, with the one of interest 
located at the distance where we see the CMB. This can 
be calculated using standard FLRW codes, but with the 
line-of-sight parts corrected for [301 EH HE] • 

For line-of-sight effects, we need to use the full void 
model. These come in two forms. The simplest effect 
is via the background dynamics, which affects the area 
distance to the CMB, somewhat similar to a simple dark 
energy model. This is the important effect for the small- 
scale CMB. The more complicated effect is on the largest 
scales through the Integrated Sachs- Wolfe effect (see [46] 
for the general formulas in LTB), which we don't con- 
sider here, since it depends on the detailed evolution of 
perturbations during the curvature era which are not yet 
understood in the context of inhomogeneous models (al- 
though some progress has been made [M] |3H [SS] ) . 

The CMB has been considered in void models be- 
fore using these approximations. Several papers have 
computed how the position of the peaks constrain the 
model [H[IFJ[331[SJ]. In g2J, the full spectrum on small 
scales (say £ > 100) was shown to fit essentially per- 
fectly, following from earlier work by [40] . These papers 
found that, although the CMB can be accommodated, 
it required either a very low Hubble rate at the centre 
or very high curvature for the simplest asymptotically 
flat void models, so ruling them out. Models which were 
not asymptotically flat also had problems [32]. However, 
in [52], it was shown that a varying bang time could be 
introduced to make the CMB fit for sensible central val- 
ues of measured parameters. 

More recent work 60, 61 , released simultaneously with 



the first version of this paper, used these results and com- 
bined them with other data sets to constrain the mod- 
els [BU] or even rule out void models altogether [6T| (when 
the bang time is homogeneous). The main constraint 
comes from the need of a small value of Hq for compati- 
bility with the CMB. These analyses both used methods 
for fitting the CMB which, we shall argue, over-constrain 
the models by effectively fine-tuning the radiation profile 
unnecessarily. 

A common theme in these papers was that the influ- 
ence of radiation was largely neglected and treated as 
a test field in an LTB void model when performing the 
matching to an early-time FLRW solution. This seems 
reasonable because radiation only contributes about 10% 
of the energy density at last scattering. Furthermore, the 
area distance-redshift relation only varies by at most a 
few percent in FLRW models at z ~ 1 100, depending 
on whether radiation is included or not. It seems sensi- 
ble that any errors induced by this approximation would 
lead to percent-level errors in any subsequent parameter 
estimation. 

We shall argue here that this is not the case. By in- 
cluding the effects of radiation in the dynamics of the 
spacetime, we find that the CMB constraints are con- 
siderably less restrictive, even for a high value of the 
Hubble rate at the centre (i.e., h ~ 0.7). We also find 
that asymptotically open models are preferred, not closed 
ones. The effects of radiation are important, as first sug- 
gested in [58], an analysis we expand upon in detail here. 
There are several (interconnected) ways to see this: one 
is simply that the extra degree of freedom from including 
inhomogeneous radiation absorbs a constraint which, we 
argue, is artificial in pure dust LTB (as we describe in 
more details at the end of Sec. Ill]); another is that the 
dynamics of radiation at early times imprints itself in the 
dynamics at late times (see Appendix A 1 1 ; a third is that 
CMB observations cannot determine, on their own, the 
interior of our past lightcone, which appears as a reper- 
cussion in other analyses. 

In previous analyses an LTB relation was used in go- 
ing from scale factor to redshift, a function a{z). In LTB, 
this differs from the FLRW relation a(z) = 1/(1 + z) by 
tipically a few percent out to z ~ 1 100. A direct con- 
sequence of this relation is that the CMB temperature 
today, on a surface of constant time, varies by at most a 
few percent across the void for a monotonic profile. Con- 
sequently, the variation in the radiation density must be 
homogeneous at the sub-percent level. In effect, this is 
the condition which requires a very low Hubble rate at 
the centre. Does this make sense? Intuitively it does 
not because the difference in the Hubble rate inside and 
outside the void at late times varies typically by tens of 
percent, and the density by even more than this - fac- 
tors of 5 or more are not unreasonable. So, inside the 
void, an observer is expanding away from the CMB sig- 
nificantly faster than an observer outside. Surely they 
would measure a significantly cooler CMB than an equiv- 
alent observer outside - by much more than the percent 



level? 

More concretely, consider a void much larger than the 
Hubble scale today, so that we can safely use FLRW re- 
sults inside the void, and asymptotically far from the 
void. If the model has homogeneous baryon fraction and 
baryon-to-photon ratio (say, as an example to set the ra- 
diation to matter ratio) , then the CMB will be emitted at 
constant temperature, and constant Hubble rate. Today, 
the temperature in any region where we can neglect the 
effect of inhomogeneity must locally satisfy Tq oc Vt m h 2 \ 
depending on the density of the void, this can vary by a 
large fraction from inside to outside. 

As a final way to see the freedom inherent in the radi- 
ation, consider that we must be able to specify a spher- 
ically symmetric model as a Cauchy problem in the fol- 
lowing manner. We have the freedom to specify both the 
radiation density and matter density as independent ar- 
bitrary profiles on a hypersurface of constant time today, 
and evolve them backwards to the last scattering surface. 
The Hubble rate today is also a free function, but this can 
be specified iteratively to make the bang time homoge- 
neous if desired. This is a perfectly feasible way to make 
a model (in principle), and demonstrates that the tem- 
perature today is not forced to be near-homogeneous - it 
can be anything, and is independent of the matter profile. 
As far as we are concerned, if these profiles are chosen 
such that the temperature of the CMB at the centre is 
as measured (with all other central parameters chosen 
as desired), and parameters asymptotically have the cor- 
rect baryon fraction and baryon-photon ratio (which fixes 
the CMB peaks), then we have a candidate model for 
the CMB; as we discuss, only the area distance provides 
an additional constraint (to place the peaks in the right 
place). Once these are fixed the model can be evolved 
backwards to find out the interior of past lightcone of 
the central observer - it is not fixed a priori. This argu- 
ment is really just stating the obvious: that observations 
of the CMB cannot, on their own, tell us about conditions 
in the interior of our past lightcone. 2 

To use the small-scale CMB to constrain void models 
in full generality actually requires a two-fluid radiation- 
plus-matter solution of the field equations in which the 
two fluid are non-comoving in the radial direction (in the 
background). (Strictly speaking, we need a three-fluid 
description because the dark matter and baryons decou- 
ple well before last scattering.) As we shall discuss, this 
is quite a challenging problem in general. However, the 
physics of decoupling (which we assume can be described 
as in FLRW) unambiguously fixes the baryon fraction 
and baryon-photon ratio at last scattering. Then com- 



2 There is an important exception to this statement given by per- 
turbations in the observed CMB spectrum induced by microwave 
photons rescattered in clusters or inter-cluster medium which can 
in principle (together with baryon acoustic oscillation measure- 
ments) probe the inhomogeneity profile in the interior of our past 
lightcone, as we will discuss in Sec. IV. 



bining these constraints with bounds on the area distance 
(which encodes line-of-sight effects), the late time matter 
void profile and the radiation profile are constrained. For 
the latter, we shall argue, by analyzing the field equations 
using the approximation that the radial velocity between 
the matter and the radiation is small, that models can 
be constructed which can fit the CMB, alongside other 
observations (SNIa and Ho). However, given the uncer- 
tainty in the theoretical description, a full likelihood anal- 
ysis is premature, and not considered here. 

This paper is organized as follows. First we review the 
CMB in FLRW, extracting the important parts for our 
analysis. This will show how to disentangle effects asso- 
ciated to physics of decoupling from the evolution history 
of the Universe in the small-scale CMB we observe. We 
then provide a detailed prescription for calculating the 
small-scale CMB in void models (the first peak and be- 
yond) , developing the outline first presented in 58J . Con- 
straints on void models are presented in Sec. |IH[ which 
is followed by a discussion on implications for other ob- 
servations. Sec. |V| concludes. We dedicate the Appendix 
to describe the importance of including radiation in the 
analysis, which involves a full two fluid solution to the 
field equations. 



II. THE SMALL SCALE CMB - A MINIMAL 
PRESCRIPTION 

A. The FLRW case 

Before we begin to discuss a void model in which to 
compute the CMB, let us first describe our approach, 
and re-consider the FLRW case. In a spherically sym- 
metric universe, an observer at the centre observes the 
CMB at a uniform area distance (also called angular 
diameter distance) cIa(z*) in all directions, at redshift 
z*. Roughly speaking, the temperature fluctuations ob- 
served have two contributions: local contributions from 
the patch of the last scattering surface observed, which 
depend purely on the physics of decoupling (plus initial 
conditions set by the primordial power spectrum), and 
non-local line of sight effects. Let us examine some use- 
ful quantities for describing the CMB in FLRW, trying 
to separate the local from non-local effects, and, in par- 
ticular, scale factors from redshift. 

The CMB shift parameters 
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are sufficient to characterize the key features of the first 
three peaks of the CMB [531 |M] (see also [55 - 68 who 
use different but analogous parameters). Given standard 



thermal history and matter content, the physics which 
determine the first three peaks also fix the details of the 
damping tail [M] • With the exception of d,A (z*) , all quan- 
tities are local to the last scattering surface of the CMB 
that we observe, and the evolution of the universe before 
decoupling. 

Let us consider l a first. The comoving sound horizon, 
r s (a*), is given by 
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Here, the Hubble parameter is given by 
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the scale factor a is defined as unity today, and all 
f2's represent today's values; we also define Hq = 
100/ikms~ 1 Mpc~ 1 . The full radiation contribution is 
given by 
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where N e g is the effective number of relativistic degrees 
of freedom around today - the standard model has N c g = 
3.04. 

Now, a* = a*(t*) is the scale factor at the time of 
decoupling, t*. Hence, the proper size of the sound hori- 
zon is a^r s (a if ). This means that l a is the ratio of area 
distance to proper distance, which are both well defined 
distance measures in any spacetime - unlike comoving 
coordinate distance which has no physical meaning in 
general spacetimes (for example, WMAP jSHKHT] use co- 
moving distances rather than physical distances). 

Now consider l oq . The equality scale is given locally 
by the time of matter-radiation equality, when a oq = 
£l r /£l m . The mode which enters the Hubble radius at 
equality has comoving wavenumber 
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The proper size of the Hubble radius at equality is just 



H cql so it's useful to think of l cq as 
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again, a ratio of area distance to proper distance. 

For variables which are local to the last scattering sur- 
face (LSS) it is possible to write them in terms only of 
local quantities. In particular, we can remove all ref- 
erence to 'today' where a = 1 and the fi's are defined 
(or whichever reference time is chosen). Changing vari- 
ables from scale factor to temperature using a = T /T, 
where To is today's temperature, we find that all refer- 
ence to H , T and the fi's, etc., factor out of the quanti- 
ties above, leaving just the baryon fraction /& = Q(,/f2 m , 



the local temperature of decoupling T* , and the baryon- 
photon ratio r\ which defined as the number of baryons 
per photon, 
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Thus, the proper size of the sound horizon at decou- 
pling may be written as 
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where H/Hq f=a 2998 Mpc and we have defined the 
dimension-full constants 
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Crucially, these have no dependence on any parameters 
of the model. That is, we are not free to specify the 
tn's, apart from N c g. These are derived assuming that 
jb and r\ are constant. 3 Hence, the proper size of the 
last scattering sound horizon, a*r s (a*), depends only on 
i], fb and T*. The temperature of decoupling at the level 
of approximation we use also only depends on n and /&: 
The elastic Thompson scattering rate depends only on 77, 



r(T) 



\/gC(3W 1/2 /mg \ 3 / 4 T 3 -A/2T (m 

■k ' \2ttTJ ' V ; 



which at decoupling must equal the Hubble rate 

iook„^Lpc-. -v /( °'- + ^ )r,+ ° t i T3 ' (13) 

which also only has dependence on the local parameters 
r\ and /;,. (Using the approximation that r(T, t ) = H(T*) 
overestimates the temperature by 2 — 3% - we use the 



3 A similar expression was given in I68| using redshift instead of 
temperature. Although we don't have dimensionless constants 
doing it our way, temperature along the worldline up to last 
scattering is well defined by the local radiation density; redshift is 
ambiguous in an inhomogeneous universe as it applies along null 
cones, not timelike worldlines. Furthermore, redshift requires 
a fixed observation point which we shall find useful to remove 
where we can. 



Saha approximation for illustration only; see below for 
the more accurate method we actually use.) For l eq , the 
scale which is compared with the area distance to the 
CMB is 



Finally, 
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Thus, apart from <1a(z*), the key features of the CMB 
power spectrum constrain directly /{,, 77 and T* along the 
radial worldline of the patch of the LSS we observe. Given 
an observation temperature of the CMB, these translate 
into £l m .b,rh 2 using Eqs. (19j) - (11). We have taken the 
time to make this explicit Decause the observation tem- 
perature is not obvious in void models: that is, we shall 
find that the temperature entering in those equations is 
not the temperature observed at the centre of the void, 
but rather the temperature that a fictitious observer far 
outside the void would measure. 

Instead of the Saha approximation, we will use the fit- 
ting formula of [70] as a more accurate means to calculate 
T*, adjusted for our purposes. This is normally given as 
a formula for z* in terms of fl^h 2 and fl c h 2 , and derived 



for To = 2.725 K. Using Eq. ( 11 ) we can recast it in terms 
of r\ and /&: 
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where 
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and we have defined 7710 = 10 10 ?7. 

The parameters l a ,l cq ,R* capture all the information 
we need from the small angle CMB (as already men- 
tioned, we are not in a position to evaluate the large- 
scale features, which includes the late time Integrated 
Sachs Wolfe effect and so full calculation of the growth of 
structure). Other 'orthogonal' parameter sets are some- 
times used. For example, WMAP and others [65-68 use 
l n , R and z* , where 
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So, this relies on T , which Z Gq does not. For our purposes 
this is not good: we shall find that Tq varies considerably 



on a surface of constant time in void models, which would 
make this shift parameter awkward to use. Also, the 
other parameter in this set, z*, is not as neat for us as 
T* because it requires a point of observation, rather than 
just local physics. 



B. Inhomogeneous spherically symmetric universe 

As we observe from the above equations, the local part 
of the CMB constrains only local variables which are 
meaningful at the time up to decoupling: 77, /&, and T* 
which is fixed by 77 and fb (or, alternatively, one can ob- 
viously choose any orthogonal combination of such vari- 
ables). Provided the FLRW approximation is valid be- 
fore decoupling in a horizon size then the equations above 
which involve ro's rather than Q's are sufficient. There 
is no reference to late times in the equations and it does 
not matter if the spacetime is FLRW after decoupling. 

Let us fix fb and 77 at the radius of decoupling such 
that i?* and l a /l e q satisfy observational constraints. This 
implies T* from i?* or better from Eq. ( 16 ). Hence, H* = 
H(T*) is given unambiguously by Eq. (13). In addition, 



if an observer measures the CMB to have temperature 
T obs , then the redshift of the CMB will be 1 + z* = 
T*/Tg bs , since the black-body spectrum is conserved in 
any space-time [7T]. The final constraint from the CMB 
parameters then must be d^(z*), which follows from l a or 
l C q. For a spherically symmetric universe observed from 
the centre, d^z*) is just one constraint on the model. 
The current WMAP7 constraints on the three degrees of 
freedom are illustrated in Fig. [T] A slight degeneracy is 
present between /& and dA- Observational data require 
7710 « 6.2, f b « 0.17, and d A ~ 13 Mpc (which will 
be used as benchmark values in some of the following 
figures). 

We now want to link such constraints to the param- 
eters of the model. For our purposes, we can calculate 
nearly everything along the central worldline and along 
one at the radius of the CMB (r 3> void width) , where 
we assume an FLRW description holds with separate pa- 
rameter values at the centre and asymptotically, denoted 
(in) and (out), illustrated in Fig. [2] The exception is the 
area distance to the CMB, as discussed below. 

The parameter values labelled '(out)' are, strictly 
speaking, asymptotic values and are not necessarily the 
values of the parameters along the worldline where we 
observe the LSS (which we refer to as the 'CMB world- 
line' located at coordinate distance r„). Since we consider 
inhomogeneities which are Gpc in size, and the LSS is lo- 
cated where the profile is flat 4 , decoupling occurs in the 



4 There may be effects from the slight radial gradient at decou- 
pling. The profile can be chosen exactly flat at r* if required. 
As we are neglecting the bang time function, any other type of 
inhomogeneity should also be small at this time. 
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FIG. 1: Constraints on /;,, r\ and (1a(z*) derived from 
WMAP7 constraints on a ACDM model. The points are 
sampled assuming Gaussian errors on l a — 302.44 ± 0.8, 
l cq = 137.5±4.3, WOflth 2 = 2.258±0.057, z* = 1090.79±0.92, 
the last two of which are used to calculate R* . The variables 
are computed using the fitting formula for T* , Eq. {16) , in- 
dicated by the colour. Note the slight degeneracy between 
fb and (Ia{z*), while the other combinations are essentially 
independent. 



asymptotic FLRW region. However, an observer at that 
radius will eventually see the inhomogeneity, and so the 
FLRW approximation will no longer be valid. 

Our approach relies in some respects on the spacctime 
at the centre being locally (very close to) FLRW, but 
we do not rely on this for our key results. We give a 
discussion of this in the appendix. 



C. The CMB calculation for a void 

Voids may be described by a profile for the present-day 
matter density parameter. We shall typically consider a 
Gaussian as a benchamrk case (although the void shape is 
only mildly constrained by SNIa at low redshift, e.g. \$*" 



and we discuss its impact on obtained result in Sec. Ill B ) 
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and we can assume a similar profile for the radiation den- 
sity Q r (r) (usually ignored but important here). The 
Hubble rate is given by a generalized Friedmann equa- 
tion, which we discuss in Appendix [XI When in the mat- 



ter era and radiation can be safely ignored this is given 
by the usual LTB relation (see, e.g., [30]). 

Along the central worldline and along one at the radius 
of the CMB (r ^$> a), we can neglect all r-dependence, 
and use the standard Friedmann equation, Eq. ([3]) with 
separate (in) and (out) parameter values. The area dis- 
tance to the CMB is instead calculated in the normal 
way for an LTB model 30, 48 , and is affected by the 
shape and width of the void. The CMB places constraints 
on the void model through Eqs. ffity - (111, which hold 
provided the model is asymptotically FLRW, as well as 
the area distance to the LSS. The remaining part of the 
model is specified by local observations of the CMB tem- 
perature, the Hubble rate, SNIa observations and so on. 

To fully fix the model we need to specify the parame- 
ters 



To, h, Q m , f b , r), N eS , 



(22) 



at the centre and asymptotically. We choose a(to,r) = 1 
to fix our gauge. We constrain the void as follows: 

Inside: the centre of the void 

Core constraints: 

• From the CMB temperature today we have 



T, 



(in) 



J S ""(in) 



2.469 x 10" 



)(in) 



1.690V for iV eff = 3.04 from 



2.725if = 

This gives Q, r m ' 
Eq. ft}. 

• The local expansion rate should be taken from local 
observations only, unless a full likelihood analysis 
is performed. For example, [7J] find h^ in "> = 0.742± 
0.036 from low redshift SNIa. 

• A void model which fits the SNIa well can have very 
low density at the centre. We can expect 

n^ ~ o.i - 0.2 

from SNIa constraints. 5 
Additional parameters: 

• The baryon-photon ratio can vary radially, reflect- 
ing the initial radiation distribution. From 7 Li 
measurements we can set (for models such that r\ 
is conserved or nearly conserved) [58] 



V 



(in) 



5x 10" 



alternatively, we can fix 77 as a constant and de- 
termine it directly from the CMB. Note that r]( m > 
refers to the baryon-photon ratio at late times at 
the centre, which is not the primordial value if r\ is 
not conserved (see the Appendix). 



5 Age constraints favour a slightly lower value of h than the 1721 
value, but with f2™ this low, recent constraints 1571 are rela- 



tively easy to satisfy. We don't consider this further here. 
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FIG. 2: Schematic of our model and assumptions. The central observer sees a small patch of the last scattering surface, which 
we assume evolves as in the asymptotic part of the model which is FLRW. The shaded regions are not required to be fixed for 
the calculation of the small scale CMB, but are constrained a posteriori. 



From Eq. ( 11 ) we can find 



n 



(in) _ 



VJ h - 



7 (in) T 



hi V 

(in) 



(in)* 



This fixes the local baryon fraction at the cen- 
tre / b = $% n /Om ■ An obvious choice is to 
make the assumption that /t, =const. so that 
/ & (in) = / 6 (out) . Given that the CMB constrains 

f b , this puts an additional constraint on ilm ■ 
Local observations of fa could instead be folded in 
at this stage; for example, from observations of gas 
in clusters one finds fb ^ 0.11 [73], with a spatial 
variation < 8% up to z ~ 1 [74] , 

These last two are not necessary for the CMB analysis, 
but are useful to define nonetheless. 

Now to make the link with the asymptotic parameter 
values, we calculate 



t 



(in) 



da 
~a~H 



along the central worldline. We assume the Bang time to 

be homogeneous is = 0, so ig — *o ^ s ^ ue a § e °^ ^ e 
universe everywhere today (the bang time function is a 
decaying mode in LTB models [75]). This fixes the time, 
at a(to) = 1, where the asymptotic values of quantities 
are defined. 

Note that the only parts of the calculation where we 
have used the assumption that the spacetime along the 
central worldline is FLRW lies in calculating the baryon 
fraction on the inside, A , and the overall age of the 
model. Thus, in cases where this is not satisfied, this 



will require an adjustment of /> 
give below. 



(in) 



from the values we 



\^) Outside: asymptotic parameter values 



Given values of l a /leq arL d R* we find values for 

^(Out) and ^(out)^ and; hence ^ rp(oUt) from Eq ^ 



The redshift of the LSS is 



1 + z* 



T t 



(out) 



,(in) 



(24) 



If we ignore the contribution from radiation, the 
age of the universe, t Q ou 



fixes ft 



(out) 



given /i(° ut ) 

For arguments sake, let's 

(out) and fix M out ) from this. 



or /j(° ut ) given il (out) 
assume we choose fl 
(The radiation density makes virtually no differ- 
ence to the age.) 



The observed area distance to the LSS, d^ s (z*), 
derived from l a or Z cq , may be approximately cal- 
culated in an LTB void model, neglecting the effect 
of the radiation (we discuss the errors from this ap- 
proximation in the Appendix) . This constrains the 

void profile as well as flm (given /i( out ) from to) 
but fixes neither uniquely 

Correlation with SNIa data may be used to break 



the degeneracy on a and £2 
FWHM ~ 6Gpc and fl m (z 

a£ n) ~o.i @g.) 



(out) 



(SNIa favour 
1) ~ 0.5 when 



We now calculate Tq OU ' from Eq. (|ll| . We can now 



*(out 



find Q\, r trivially (given N c g). Recalculating the 



age can now refine the estimate of /j( out ) . (We find 
these precisely iteratively.) 

A final constraint arises from the measurement of T 



-j(out) 



and its link to T* ; through Eq. ( 24 ) , namely when we 



integrate out along the past lightcone from the centre out 
to z* . In terms of time rather than temperature, it says 
that the local time at that point must equal the time 
obtained by integrating up along the timelike worldlinc 
from the big bang up to decoupling. That is, 



to — t*{z*) 



dz 
1 + z 



(d t lnx) 



nullconc 



(25) 



where x(*) r ) = —dt/dr evaluated on the past nullconc, 
and £*(z*) is the local time of decoupling at the redshift 
observed from the centre, which must be equal to 






dT 1 
Tr ~T H{T) 



(26) 



where H(T) is given locally by Eq. ( 13 ) (with no ref- 
erence to late times). Because the integral in Eq. (25) 



traverses the void, it is a constraint on the profiles. This 
constraint is extremely sensitive to the details of the ra- 
diation inhomogeneity (and the 'out' parameter values) 
because most of the integral cancels the to on the l.h.s.. 6 
Indeed, since t*/to ~ 10~ 5 , we need to know the whole 
spacetime to better precision than this; matching an LTB 
void to a radiation filled FLRW model cannot achieve this 
accuracy (see the Appendix). If the central worldline is 
FLRW, then this constraint must be automatically satis- 
fied and can be inverted to find the density profile of the 
radiation between (in) and (out). For the models consid- 
ered here, this typically gives an inhomogeneity 0(1) in 
the radiation density at early times (see Fig. 2 of [58]). 
We expect that as the void size decreases this will lead 
to stringent constraints on the radiation inhomogeneity; 
but for the Gpc scale voids we are interested in we expect 
the constraints to be mild since the void scale is only en- 
tering the Hubble scale today. Essentially, this will tell us 
how to map the initial radiation profile onto the central 
worldline evolution for the temperature. 



However, Eq. (25 1 deserves a more detailed discussion, 



which is presented in Appendix [A] Indeed, one of the 
main differences between this work and previous analy- 
ses resides in the introduction of inhomogeneous radia- 
tion which can significantly affect this £*(z*) relation (or, 
equivalently, a*(z*)) with respect to the relation com- 
puted within the LTB framework matched to a radiation 
filled FLRW model. In essence, z* is fixed by the CMB, 
but £* and a* depend on whether radiation is included 
or not. This changes things considerably. Indeed, con- 
sidering Eq. ( 25 ) as a constraint on to given z* and t* , 



Indeed, if we evaluate this constraint in FLRW, t„ varies by tens 
of percent depending on the amount of radiation present. 



this effect has similar consequences to the introduction 
of an inhomogeneous bang-time function since t* is the 
time since bang in that context. 

All the key features of the void are now known and we 
can derive quantitative constraints. 



III. RESULTS 

A. Constraints from physics of decoupling 

First, we consider constraints arising the from physics 
of decoupling, namely from l a /leq and i?* parameters. 
Here we do not take into account the bounds from the 
area distance. As an illustration of how the CMB di- 
rectly constrains the model, consider the left panel of 
Fig. [3] where we take an asymptotically flat void with 
/j( m ) — 0.7. Assuming both r\ and fb constant, we see 
that the CMB constraints translate into constraints on 



For this kind of void, 



n& n) and h<- out \ as well as T (out) . 

£lm turns out to be quite large, and not easily compat- 
ible with SNIa data (moreover, the area distance is too 
low unless the void profile has a dip of very low density in 
it, which would be again tough to fit to SNIa; see below). 
Releasing the assumption of a constant baryon fraction 
or constant baryon-to-photon ratio, one can fix the mat- 
ter density inside. In the first case (middle panel, now in 
a model which is open asymptotically), the derived /^ 

is significantly larger than f^ and in possible tension 
with cluster data [74]. When 77 is allowed to vary in- 
stead (right panel), ?yj ra 3.5 fits the CMB, which may 
also give an explanation to the 7 Li problem of standard 
cosmology as proposed in [SS] (actually, 7 Li constraints 
favour a slightly higher value than this). 

We gain a more systematic understanding of the con- 
straints implied while ignoring the area distance in 
Figs. [4l — [5| We choose some CMB parameters l a /l e q,R*, 

which imply / fc and rj^ out \ and show how these con- 
straints translate into the parameters of void models. 

Consider the case when fb and r\ are both constant, 
as in Fig. |4| As we vary h^ n \ flm is fixed by the as- 



sumption on fb through Eq. ( 23 1 (top panel) . Other pa- 



rameters such as Tq ' and ft,( out ) depend on flm ' and 
ft( m ). Indeed, /i( out ) is set by the age t ° u (assuming a 
homogeneous bang time). The bottom panel shows that 
larger M ln ) means shorter age and so larger h^ out \ and 

that larger curvature ^Jf" (i.e., lower fim ) would in- 
crease the age and so needs to be compensated by a larger 
^(out) y\iq Tn dependencies (central panel) follow 



again from Eq. ( 23 1 , namely, radiation density scales lin- 



early with matter density when fb and 77 are constant. 
In Fig. [5] we release the assumption that 77 ^constant. 



Increasing r/ m ) implies, for a fixed h^ and T , a larger 
fim ) which leads to a shorter age, and, in turn, for a 
given Hm , to higher /i( out ) (lower panel) and T 
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FIG. 3: Constraints on an asymptotically flat void, ignoring the area distance to the CMB (left). Models which fit l a /leq and 
R* are picked as in Fig. dTp, which constrains /j and r), which we assume spatially constant in this case. Because r\ is constant 
rather than fixed at the centre, this results in a scatter in the void parameters shown. In the middle, we have chosen an 
asymptotically open void, which lowers Tq , and we have released the assumption on the baryon fraction; instead we have 
chosen the matter density inside as fixed. The baryon fraction at the centre now takes up the scatter in values from the CMB, 
resulting in a fixed void with fixed h^ out ' (derived). Similarly, we can choose the baryon fraction fixed and let ry> iu ' take up the 
scatter - shown right - in which case these parameters favour a low value of ?/' , in line with the lithium constraints |58| . 



h (c 



■j (out 



(top panel, see Eq. (23) for fixed / b (out) and 7?( out )). The 



vs. flm relation is given by the age constraint 



as in Fig. Ul and T vs. flm is again related to 
Eq. pi. 



B. Constraints including the area distance 

Now let us consider the area distance to the CMB. 
This mainly depends on the void profile and the Hub- 
ble rate at the centre (and z*). In Fig. [6] we see how, 
changing ft/ ln ), flm , fim and the width a. In essence, 
the lower the density, the higher the curvature, and the 
larger the area distance to a fixed redshift. The longer 
the null geodesic spends in a low density region, then, 
the further the CMB will be located. The CMB con- 
straint of (Ia(z*) — 13Mpc implies (with the caveat that 
we are assuming cIa not to be affected by radiation inho- 
mogeneity) that for a high Hubble rate at the centre the 
matter density asymptotically must be less than unity - 
an asymptotically flat void is very difficult to achieve, in 
agreement with [40l |42] - and the width must be several 
Gpc (also implied by SNIa). Lowering the central Hub- 
ble rate pushes up the curves in Fig. [6] For a very low 
central Hubble rate a flat or closed model asymptotically 
is required to fit cIa{z*)- 

To summarize, then, let us consider how to use these 



figures to estimate a specific case. Let's choose h^ — 
0.7, and t)\_ = 4 to coincide approximately with lithium 
constraints. Then, if the baryon fraction is constant, we 
observe from the lower panel of Fig. [5] that we must be 



on the curve with f2 



(in) 



0.17 (which is fine for SNIa). 



Turning to Fig [6) we can see that for this matter density, 
and, say, a FWHM ~ 2Gpc, we must have f2 



(out) 



0.7. 



(out) 



Returning to Fig. R we have /i( out ) sw 0.58 and T c 

3.3. A void with these parameters will give a decent fit 

to the CMB. 



In the recent work of Refs. [60|l61j. it has been found 
that the CMB constraints are satisfied for asymptotically 
closed models with h^ < 0.5. It's not difficult to ex- 
trapolate this case from Fig. [6] and to note that we also 
find the area distance is correct. Although such models 
are not specifically included in Figs.|4]-[5l they are con- 
structed such that / 6 w 0.17 and 7]{ 3 u w 6.2, which 
implies, given what we said above, that they also satisfy 
constraints from physics of decoupling. Therefore, fol- 
lowing our approach, these models constitute a subclass 
(i.e., with nearly homogeneous radiation, see Appendix) 
among the void models which can fit the CMB (and we 
do not particularly emphasize them mainly because of 
the low ft/ m ) which is tension with data {72]). 
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FIG. 4: Constraints on a general void as a function of the 
Hubble rate at the centre, ignoring the area distance to the 
CMB. Here we have chosen the CMB parameters as indicated 
which imply fb — 0.17 and T)io = 6.2, and we have chosen 



these variables as constant. Different choices of Q. 



(out) 



then 



give different h^ out ' and T , but ilm (top) is constrained 
purely by the assumption on the baryon fraction. Note that 
models with T (out) ~ T (in) have very low ft (in) . 



An asymptotically flat void? 

A key finding of previous work is that the CMB rules 
out asymptotically flat voids [301 H2]) unless the bang 
time is varying, or the matter density or Hubble rate are 
much lower than otherwise observed. Under our approxi- 
mations for the computation of d^ , we also find the same 
result, principally because an asymptotically flat model 



T, 



(out) 



fc (out) 
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FIG. 5: Constraints on a general void as a function of asymp- 
totic matter density, ignoring the area distance to the CMB. 
Here we have chosen the CMB parameters as indicated which 
imply fb — 0.17 and rf^ — 6.2. We have allowed n^ n ' to 
take on different values, which, when the baryon fraction is 
constant, implies differing Qm (indicated). 



with sensible h^ and flm has the CMB far too close 
to us, if the CMB worldline is in the asymptotic region - 
this can easily be inferred from Fig. [6j Is there any way 
to make an asymptotically flat model fit? 

It has been shown in [43J HI] that one can add features 
to the radial profile to alter the area distance to the CMB. 
For example, since SNIa probe the distance modulus only 
out to z ~ 2, there are no distance observations between 
them and the CMB. This freedom can be utilised to put 
an under-dense shell around us, thereby increasing the 
distance to the CMB by the required degree. 

There is another way to make asymptotically flat mod- 
els fit. If we increase the number of relativistic degrees 
of freedom above the canonical value of N c g = 3.04 
(and allow for a spatially varying /;,), then constraints 
in Fig. [l] will be modified and, in particular, this can 
bring down the area distance to the CMB. We find 
that with N g ~ 8 an asymptotically flat model, with 



h( 



0.7,0 



(■") 



0.2,FWHM w 4Gpc fits the CMB 
perfectly well. This is a similar conclusion to |39j and 
deserves further investigation in this context. 

These scenarios are depicted in Fig. [7] 
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FIG. 6: Area distance to the CMB as a function of asymptotic 
matter density. Here we have used CMB parameters such that 
z» = 1090. Clearly, there is no problem making a void model 
with the correct area distance to the CMB, even with a high 
value of hr, provided the model is open at the CMB radius. 
The CMB constraint is shown for N*g = 3.04; increasing this 
lowers the required area distance to the CMB. 



Pi 



cos 






Pi (cos ( 



d-AJz* 
<1a{z* 



c [d A (B,)/d A (z,)]e 



for £> 1,(27) 



where it is assumed that Cg is smoothed appropriately 
for non-integer £. 

As the final step to computing the CMB power spec- 
trum, we can use CAMB [75] as follows. First compute 
a void model which is going to give appropriate CMB 
parameters, according to the steps given above. Then 
compute the CMB power spectrum using the 'out' pa- 
rameters for h, O m /i 2 , Qbh 2 and T in CAMB (using a 
zero A model). Compute rf^C-S*) f° r this model. The 
shift parameter S = dA(z*)/dA(z*) may be calculated 
with dA{z*) found from l a (since a can be tweaked to 
make this precise). Then, the output of CAMB, which 
gives £(£+l)Cg, requires the £ axis compressed by S with 
no shift to the amplitude, since, using S = £/£, 

£(£ + \)C t *■ W + l)S 2 } [S- 2 C e/s ] ~ £(l + l)d t . (28) 

As an example, let us compare with a flat ACDM 
model which fits the WMAP7 data well. Choose l a = 



C. CMB power spectrum 

Once the void model is fully specified the power- 
spectrum of the CMB may be calculated properly using 
the procedure given in [33] [55], for large £. This is an im- 
portant final step because the observational constraints 
on the 3 parameters we have used to constrain the model, 
l a , l C q, R*, are actually derived parameters from the data, 
using an FLRW model [57] . So, a void model which gives 
particular values for l a , Z eq , i?* may have a slightly differ- 
ent CMB power spectrum from an FLRW dark energy 
model with the same l a ,l e q,R*- We find this difference 
of order a percent or so, in broad agreement with using 
the parameters for different dark energy models |65j . 

To calculate the CMB angular power spectrum for a 
given void model, first calculate the Cg's for an FLRW 
model using the asymptotic parameter values. This 
power spectrum will have the same intrinsic CMB fea- 
tures as the void model, since they see the same LSS, but 
will need to be shifted to account for the incorrect area 
distance used in the calculation. Let the area distance in 
the FLRW model be dU(£*), where the LSS is at redshift 
z». The relation to convert to the power spectrum the 
void observer will see is [68] (see also [42|l69]) 
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FIG. 7: Examples of inhomogeneous models with a central 
void providing an area distance to LSS able to fit CMB 
data. The computation of cLa has been performed assum- 
ing LTB framework. Open models (black dotted) naturally 
fit. Asymptotically flat models either requires a profile with 
a shape more complicated than a simple Gaussian void (red 
dashed and green solid) or an JV e a significantly larger than 3 
(blue dashed-dotted). 
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FIG. 8: The TT (middle) and EE (bottom) angular power 
spectra for a flat ACDM model and a void model which give 
the same l a ,l B q, and R* parameters. The difference between 
the two is just a few percent (top). Although the void is de- 
rived with the assumptions indicated, we find similar plots 
for different types of void which fit the CMB parameters in- 
dicated. 



302, l cq = 136.6 and R* — 0.63 which give a good 
fit. Solve Eqs. for the flat ACDM parameter val- 
ues. This gives: tt b h 2 m 0.0226, Vt m h 2 ra 0.133, and 
h = 0.713, in good agreement with the WMAP7 best fit. 
We show the TT and EE power spectra for this model 
with zero tilt in the primordial spectrum in grey in Fig. [8] 
Let us compare this with a void model which has the 
same parameters l a , l cq and i?*. As an example, let us 
choose a model with constant fb, open asymptotically 
with Qm ' — 0.7, with a high Hubble rate at the centre 
to fit HST data M in ) = 0.74, and a « 6 Gpc. Derived pa- 
rameters with these assumptions, and taking a constant 



7], are il; n 

(in 

summg rj 1 



(in) 



0.242 (or, alternatively, fi 

] = 4), fl ( ° ut) » 0.120, h(° ut 1 



a 0.156 as- 
0.636 and 



T, 



(out) 



3.15 K, the last 3 of which are used directly in 
CAMB. Using a shift parameter of S *=* 0.868, we find the 







CMB TT and EE power spectra shown in blue in Fig. [8] 
As we can see there is very little difference between the 
FLRW model and the void model. Nearly identical fig- 
ures are found for other void models which produce the 
same l a ,l C q, and R* parameters. 



IV. DISCUSSION 

The CMB in inhomogeneous universe models is clearly 
complicated (see the Appendix). Even when spherical 
symmetry is assumed the field equations are too difficult 
to give a tractable solution when radiation and matter 
are to be accounted for. The reason is that because the 
density is not homogeneous, pressure gradients in the 
radiation force the radiation and matter to have differ- 
ent rest-frames and become non-comoving after decou- 
pling. In fact, CDM decouples from the radiation and 
baryons much earlier, and so we really have to deal with 
a three-fluid solution to properly model the spacetime - 
and that is before perturbations are included! As we 
have discussed, these effects appear to be critical for a 
complete analysis of the CMB. Indeed, it is important to 
note that the void amplitude is significantly larger than 
other perturbations in the radiation era 58 , so details 
of the radiation era are critical for a full analysis. Given 
such uncertainties in the theoretical parameter estimates, 
we considered a full likelihood analysis involving different 
cosmological data to be premature. 

We have seen how the CMB constrains the baryon frac- 
tion and baryon-photon ratio severely at last scattering 
(as in FLRW), as well as the area distance to the CMB. 
Assuming that this can be computed using the LTB so- 
lution, this restricts the void profile to be an open model 
at large distances (> 13 Gpc say, in comoving distance), 
but is otherwise broadly in line with SNIa and local Hub- 
ble rate measurements. Models which reach asymptotic 
flatness at the CMB distance can fit, but require a shell 
around us of low density (for ft/ m ) > 0.6); whether this 
is at all plausible (i.e., not involving strong fine-tuning 
in the formation mechanism) remains to be seen. Al- 
ternatively, a high number of relativistic species around 
today appears to allow asymptotic flatness too (in line 
with [31]). Of course, much more freedom is possible if 
the bang time is inhomogeneous [32], and it is an im- 
portant question exactly how this further opens up the 
possibilities for void models which fit the CMB. We do 
not discuss this possibility in detail, taking a conserva- 
tive approach by assuming variations in the bang-time 
much larger than time of inhomogeneity creation to be 
unlikely 

We have argued that the final constraint from the 
CMB, Eq. (251, is more nebulous, since it is a constraint 
on the primordial radiation density profile (or, alterna- 
tively, the bang time function). We argue in the Ap- 
pendix that a radiation density profile with a similar 
FWHM to the matter and an inhomogeneity 0(1) around 
decoupling, which gives the correct temperature in and 
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out that we calculate using our FLRW approximations, 
will easily satisfy this constraint. However, given that 
the full two/three-fluid solution is lacking, this requires 
further investigation to quantify properly. 

Is it physically reasonable that small changes to the de- 
tails of the radiation era along the central worldline and 
asymptotically can lead to 0(1) changes to the model 
today? Yes. Even though the radiation free-streams af- 
ter decoupling, the imprint it leaves in the dynamics at 
early times is very important. At equality, if the Hubble 
rate at the centre and asymptotically are even a little 
bit different, then this difference in Hubble rate will typ- 
ically grow when measured on surfaces of constant time 
(unless the curvature is constant). (We discuss this in 
the Appendix (see Fig. 10).) Conversely, then, evolv- 



ing an inhomogeneous model backwards from today (if 
we knew how to do it) results in mild inhomogeneity at 
early times. 

In other words, we found that voids including a non- 
adiabatic mode, i.e., an 0(1) isocurvature mode at early 
times, can fit CMB data much more easily with respect 
to the pure adiabatic cases considered in the literature 
so far. 

A number of other issues can now be further under- 
stood, in light of our analysis, where properties of the 
CMB are used for constraining cosmological models. Let 
us consider briefly the BAO and measurements of the 
CMB anisotropics as seen by distant observers. 

Baryon Acoustic Oscillations: In principle, the 
BAO are a smoking gun test of inhomogeneous models 
because they provide a measurement of H(z) indepen- 
dently of distance measurements [32l [40J [41]. An LTB 

model which fits the SNIa data will have a different 
H(z) from a dark energy FLRW model with the same 
distance modulus (unless the bang time function is 
fine-tuned [45] )■ However, the BAO rely on compar- 
ing the size of the sound horizon at the baryon drag 
epoch with the same feature imprinted in the galaxy 
power spectrum. Irrespective of the details of structure 
formation, without which no fully reliable predictions 
can be made, this method suffers a serious drawback 
in void models in directly determining H(z). This is 
because the sound horizon at the baryon drag epoch 
need not be spatially constant. The BAO technique 
measures (see e.g., [77]). from the spherically averaged 
BAO peak positions at redshift z, r s (a,d)/Dv(z), where 
Dv(z) is the 'volume distance' oc [cIa(z) 2 /H(z)] 1 / 3 , and 
r s(0'd) is the comoving sound horizon at the baryon drag 
epoch a c i along the timelike worldline of an observer 
seen at z, which is inside our past lightcone. In FLRW 
this is just the same as measured in the CMB angular 
power spectrum (strictly speaking, it's inferred from 
it since Zd < z*). In LTB, this is an extra degree of 
freedom which must be measured separately, unless it is 
assumed. Even taking ratios of BAO measurements does 
not help because the sound horizon parts do not cancel 
as they do in FLRW. In principle, then, the BAO tell us 



only about the radial profiles of fb(r) and r](r), and not, 
unfortunately, H{z). (Of course, to study the BAO in 
detail structure formation must be properly calculated.) 

Off-lightcone CMB observations: There have been 
a number of suggestions of how to constrain the dipole 
and higher anisotropics of the CMB as seen by distant 
observers [3TJ [351 W2\ EH] • These look for spectral distor- 
tions of the CMB from anisotropic scattering [3T] [75] , or 
they use the kinematic Sunyaev-Zeldovich effect to mea- 
sure the peculiar velocity of clusters with respect to the 
CMB frame [35] [53] ■ These elegant methods rely prin- 
cipally on the large CMB dipole as seen along our past 
lightcone by observers comoving with the matter, which 
is proportional to the peculiar velocity between the CMB 
rest frame and the matter rest frame (see the Appendix 
for the definition and discussion of this). However, be- 
cause the temperature of decoupling T* has freedom as a 
function of r from the freedom in fb(r) and r](r), it can be 
adjusted to make the dipole seen by off centre observers 
small. We can see this is feasible because the dipole calcu- 
lated approximately in LTB comparing the redshift look- 
ing in, to the redshift looking out, is of order the percent 
level; the change in T* possible by fiddling mildly with 
fb and n is also of order the percent level (see [62] for 
an example). As we discuss in the Appendix, when con- 
sidering the two-fluid solution, to calculate the peculiar 
velocity we need to know the radiation profile (and a full 
integration of the field equations). In essence then, spec- 
tral distortions of the CMB and the kSZ effect in void 
models are a direct measurement of the radiation profile. 
While measurements of zero dipole for all off-centre ob- 
servers would imply unrealistic fine-tuning problems and 
so effectively rule void models out, it remains to be seen 
whether sensible profiles and models will be compatible 
with future observations of these effects. (Note that a 
systematic dipole measured would be a smoking gun for 
voids.) This is an important topic for future research, 
because if we can observationally show that the dipole, 
quadrupole and octopole of the CMB around distant ob- 
servers is small, then this is a powerful test of the FLRW 
models [T] [75] . This is also important to test the adia- 
baticity of the history of the universe. Indeed the extrac- 
tion of temperature T(z) is typically degenerate with the 
degree of isotropy of the sky at redshift z (see, e.g., [SO] 
where T(z) is extracted from SZ effect and shows strong 
correlation with the velocity between cluster and CMB 
frame which drives the kSZ effect). 

The future, then, of these 'Copernican tests' lies 
in cross-correlation of measurements, as illustrated in 
Fig. [9] Given enough observations we can reconstruct 
the conditions across the entire last scattering surface 
inside our past lightcone. 

As a final comment, we also note that constraints 
on our distance from the centre of the void from the 
CMB [15j 09] might change when a more thorough anal- 
ysis of the peculiar velocity is undertaken using the field 
equations presented in the Appendix. For example, since 
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spectral distortions of the CMB 
or the kSZ effect measured 
here is related to the 
temperature of decoupling 



FIG. 9: BAO and measurements of the CMB anisotropics 
around distant observers tell us conditions around the time of 
decoupling inside our past lightcone, and so help us to recon- 
struct the radiation profile. To use these probes to rule out 
void models they need to be cross-correlated with other ob- 
servations in different directions such that they both measure 
conditions at the same radius. 



the peculiar velocity depends on the radiation profile one 
can presumably construct models such that distant ob- 
servers see no dipole at all. In fact, this mechanism could 
be used to construct void models in which the Coperni- 
can problem is dramatically relieved (at the expense of 
a temporal coincidence problem) . Higher multipoles can 
then be used to place constraints on the distance from 
the centre, but would lead to much weaker bounds. 



V. CONCLUSION 

The ethos of modeling the universe with an inhomoge- 
neous model is different from the standard FLRW mod- 
els. Even where FLRW models contain quintessence flu- 
ids or modified gravity, flat ACDM is the simplest model 
common to all possibilities and serves as a fixed point 
from which to smoothly interpolate from. Extensions 
to the concordance model typically raise more questions 
than they answer. This means that, while ACDM is a 
good fit to the data, Occam's razor and Bayesian Infor- 
mation Criteria can be used to good effect. If it ain't 
broke don't fix it, so to speak. 

But if we consider the universe inhomogeneous on Hub- 
ble scales (and larger), then in some sense we must aim to 
reconstruct the density and curvature distributions etc., 
directly from observations, rather as if we were making 
a map. A priori assumptions such as asymptotic flatness 
or Gaussian profiles, however reasonable, are basically 
arbitrary. This is unfortunate in a sense because it gives 
us far too much freedom for model building for our fairly 
limited observations (even in the idealized case [JJ). An 
exception to this would come if we had early universe 



models for the creation of a Hubble scale void - such a 
model would presumably come with restrictions on the 
possibilities available, and so some hope of ruling them 
out. 7 

In this context, then, the CMB is not as restrictive as 
it might first appear. Indeed, we find it much less re- 
strictive than previous works which neglect the dynam- 
ical contribution from the radiation. Clearly it tells us 
that the universe must be nearly spherically symmetric, 
and, in the simplest interpretation, that we must be fairly 
close to the centre. Other than that, we have seen that 
it tells us the baryon fraction and baryon-photon ratio in 
one sphere around us, as well as the area distance to the 
CMB. This is just as in FLRW. Our main finding is that, 
given a monotonic void profile and a void which fits the 
local Hubble rate and SNIa data, a void which is open 
with flm ~ 0.7 will fit the CMB without any problems, 
as already anticipated in |58] . In particular, we find that 
the requirement in other analyses [IfJ [ST] that the CMB 
enforces a low Ho locally is an artifact of attempting to 
model the full spacetime of inhomogeneous matter and 
radiation a a separate LTB dust model embedded in radi- 
ation filled FLRW. Rather, we have argued that the full 
solution of spherically symmetric radiation plus matter 
is required before we can say what the precise constraints 
on void models of dark energy are. 
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Appendix A: The importance of the radiation 

A key part of our approach lies in including the radia- 
tion energy density. We are able to do this only because 
the spacetime is asymptotically FLRW; a comprehensive 
analysis of the CMB in a spherically symmetric model 
including radiation is a significant challenge. But why is 
it really important to include the radiation? 

In essence, it follows from the fact that 



(out) 



= T, 



(out) /„(out) 



can 

rp(in, 
1 



be significantly different from 1/(1 + z*) 
/TJ; out \ if T (out) is significantly different from T ( 



(in) 



7 This is not as far fetched as it might seem: there are a number of 
speculative ideas which can create large spherical features that 
exist today, e.g.. 
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This important effect is missed when the dynamics of ra- 
diation is ignored. Indeed, in a pure dust LTB model, 
one may calculate a* by choosing a void model, inte- 
grating out to z* using an LTB solution, and then cal- 
culating a* at that redshift. A key feature of the LTB 
solution for the types of voids considered here is that 
the function a(z) one finds is very close the the FLRW 
relation a(z) — 1/(1 + z), with a difference of a few per- 
cent. Now, since the spacetime is asymptotically FLRW, 
the temperature asymptotically today must be given by 



T, 



(out) 



T, 



(out) / (out) 



-i(out) 



/a* 



the standard FLRW relation. In 



•-(in) 



this case T and T could only differ by a few per- 
cent. But this seems to contradict Eq. dTTh , which also 
must hold asymptotically, and fixes T given /& and 
77 from the C MB, and O m /i 2 from the void shape 
we show in Sec 



III 



this gives a difference between T 



As 

(out) 



£ C q\ choose /i(° ut )(£ ° ut) ) = 10 5 ; evaluated at 12Gyr 
this model has h « 0.54 and fi m « 0.9996. Now set the 
initial conditions at the centre at the same cosmic time 
ieq 11 ; we fix the density parameters and tune h to make 
this precise. Unless the model is flat at the centre too 
this time will not be quite the central value of matter- 



While the scale of the void is much larger than the Hubble scale, 
we can calculate the evolution of the parameters of the centre of 
the void as in FLRW - since the void is Gpc, and the Hubble 
scale today is ~ 4 Gpc, we can estimate that FLRW evolution 
should be good up until, say, 0.1 — 1 Gyr. By this time, radiation 
is subdominant and makes sub-percent changes to the evolution 
history. Provided we are interested only in the overall dynamics, 
and not the radiation energy density (which may have corrections 
to FLRW evolution - see below), we can calculate the future 
evolution as in FLRW. 



radiation equality, but it will be close to it. We find 
significantly different future evolution depending on how 
set the initial conditions on the matter density and radi- 
ation density. That is, if we compare a model with ho- 
mogeneous radiation density at £ C q U we find a different 
model from choosing the matter density homogeneous by 
a significant amount at late times. This difference does 



not just decay away, as we see in Figs. 10 and 11 but 
grows to an 0(1) difference in the curvature and Hub- 
ble rate. This is because at a fixed time homogeneous 
matter vs radiation require different Hubble rates, which 
translates to different curvatures; differences in curvature 
is important. (Note that if we run the model backwards 
the homogeneous component does not stay homogeneous, 
of course, so these fine-tuned examples are just for illus- 
tration.) 



and Tq of tens of percent, not just a few. 

This makes sense: the CMB is emitted at nearly con- 
stant temperature. At the centre of the void the Hubble 
rate is much higher than it is asymptotically, so it would 
seem natural that the temperature on a surface of con- 
stant time should be lower at the centre because it is 
moving away from the observed patch of the LSS faster 
in all directions. Moreover, provided that the central and 
asymptotic worldlines are (very close to) FLRW and that 
T oc (1 + z) along null cones, a difference between Tq ' io 

and T can be simply translated into a constraint for ffj 
the radiation profile at a time around £*. Tens of per- 
cent give an 0(1) inhomogeneity in the density, which is 
comparable to the matter void. 

Another way to see this is to consider what is going 
on at early times. Could a small difference in the en- 
ergy density of the radiation at the centre compared to 
asymptotically really make much of a difference to later 
evolution? Perhaps surprisingly, yes. Let us illustrate 
this as follows. 8 Consider starting the evolution of a 
void model from matter-radiation equality to late times, 
and let us assume an asymptotically flat model for def- 
initcncss. Asymptotically at matter-radiation equality, 
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FIG. 10: Evolution of the Hubble rate for models set with 
different initial conditions at £cq" • We have chosen S — 10~ 4 
and h — 10 5 , and considered different cases: inhomogeneity 
in the matter only, inhomogeneity in the radiation only, and 
equal inhomogeneity in the radiation and matter. Clearly 
these lead to very different behavior at late times. We have set 
to — 12 Gyr. Note that we can expect this FLRW evolution 
to be accurate while the void is larger than the Hubble scale - 
that is, until the Hubble scale reaches 1 Gpc or thereabouts. 
After than, we might expect the radiation temperature to 
evolve differently, but not the overall dynamics shown here. 



Hence, radiation must be important when computing 
relations at high redshift as we now discuss. This features 
in the constraint given by Eq. ( 25 ) . 
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FIG. 11: Central curvature and Hubble rate at 12 Gyr for 
an asymptotically flat model, as a function of initial inho- 
mogeneity which is set at matter-radiation equality. The 
dashed curves have homogeneous matter content initially, and 
so the horizontal axis represents the inhomogeneity in the 
radiation; the solid curves represent the reverse situation. 
We have a significantly different curvature at late times de- 
pending on whether the matter or radiation are chosen ho- 
mogeneous. Obviously even larger differences appear for the 
generic case when both matter and radiation are inhomoge- 
neous (not shown). 



1. Using matching to understand the unexpected 
sensitivity of including radiation 

The main difference in our approach from |40j and oth- 
ers is that we relieve a constraint on Tq spatially which 
translates into a very low h^ ln K We can see schemati- 
cally why this constraint cannot be used by considering 
matching an LTB model to an FLRW with radiation at 
some intermediate redshift z m . It is assumed in [40] that 
we can choose z m high enough to be outside the void, 
but low enough for radiation to not be important. How- 
ever, matching a dust model to a radiation filled FLRW 
model typically brings in percent changes to the asymp- 
totic parameter values, as illustrated in Fig.[l2j these are 
significant. 



Consider the constraint Eq. (25), which takes us from 



the top of the central worldline today to the point of 
last scattering, as seen from the centre, along the null 
cone. This is the constraint used in other works to show 



.(in) 



(out) 



that models with radiation must have Tq 111; ~ T 
Here we will demonstrate that, without knowing the full 
spacetime including radiation everywhere, it is far too 
sensitive to the radiation era to use in any approximate 
manner. 
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FIG. 12: Matching a dust LTB model which is asymptotically 
flat, with fe (in) = 0.65, fi^ n) = 0.2 and FWHM= 3Gpc, to a 
radiation filled FLRW at some intermediate redshift z m ■ Even 
though the radiation is sub-dominant after decoupling, this 
matching incurs 0(%) level errors to the parameters from the 
asymptotic LTB model, which translates to large errors on 
the age (top). This is far too high for the accuracy required 
to determine zl /z* — T /T correctly. Matching 
is performed such that the Hubble rates agree at z m in the 
LTB model: H^ TB (z m ) = H FL , but we have to make the 
choice of exactly how to match up to the FLRW spacetime. 
We try two choices for where to match H Fh , at z = z m and 
a = aV rB (z m ); even this minor change gives large changes in 
model parameters. 



We can write Eq. (251 approximately as 

.(out) _ 



En — C* 



dz 
(1 + Z)H {] (Z) 



(Al) 



where we have assumed a simplified form for this rela- 
tion for illustration (namely, to follow other work, while 
actually this integral has additional corrections to H\\ (z) 
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from evolving curvature terms which we discuss below). 
Here, H\\ (z) represents the real Hubble rate in the space- 
time with matter plus radiation along the past null cone 
of the observer at the centre (which we don't know how 

to calculate), and z; in ' = T*/Tq — 1 is the redshift of 
decoupling as seen from the centre. The times on the left 
are the proper age of the whole spacetime (calculated in- 
cluding radiation) to, and the proper time of decoupling 
where it is observed, in the asymptotic part of the space- 
time, i* . Now, an asymptotic observer will measure 
the same time difference to last scattering but at a dif- 
ferent redshift z* = T*/Tq ' — 1, with a Hubble rate 
along their past lightcone H^° ut \z), implying 



t -t 



(out) 



dz 



(l + z) J ff(° ut )( z ) 



(A2) 



Now let us assume that we can model the matter and 
radiation model of the void as a dust model out to some 
matching redshift z m , and as an FLRW model with ra- 
diation from z m to z* , ignoring the subtleties involved 
in the matching. This implies that 



dz 
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FIG. 13: The function t(z) is insensitive at high redshift. 
Uncertainty in this integral by factors of t* result in differences 
of hundreds in redshift. Here we have shown dust models (all 
of them), and some extreme radiation + matter models for 
comparison. 



Gathering these results together, Eq. (All may now be 
written as 



The first of these integrals is the difference in time at the 
centre of an equivalent LTB dust model and the time at 
the matching point in the same model, viz.: 



dz 



(l + z)ff LTB (z) 
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The second integral may be formally split at z 
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Now, in Eq. (A2) we can formally split the integral at 



the same z m , even though it has no significance asymp- 
totically, 



dz 



(l + z)iJ(° ut )(z) 
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*o)-(4 TB -^ ut) ) = 

dz 



(°ut) (1 + Z)#(° ut )(z) 



(A7) 



The integral on the rhs will be 0(t* ) even when there is a 
difference of hundreds in the integration limits, because 
it is at such high redshift - see Fig. [13J This is a key 

result: we can expect a large difference in z* compared 
, and hence a large difference between T and 
even if the difference on the left hand side of this 

equation is 0(t*) (which is a tiny fraction of to)- in 

general we expect 



to Z: 



T, 
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{to — to) ~ t* 



(A8) 



since the age of the LTB model, which does not include 
radiation, is different from the correct age by a signif- 
icant multiple of the decoupling time along the central 
worldline (typically tens of times the decoupling time). 
Similarly, the difference in time at the matching redshift 
between the LTB dust and the matched radiation model, 
(t^ B — tm '), must also be Oft*). This implies that we 
expect, generically, 



dz 



(out) (l + z)H(° ut )(z) 



oc t* 



(A9) 
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with factors of 10 involved in the proportionality. This 
implies a difference in redshift of hundreds in the decou- 
pling redshift between central and asymptotic observers, 
which is in the tens of percent - exactly what we find 



.(in) 



-i(out) 



in Tq"' vs. Tg UUW . What this also implies is that it is 
inevitable to find differences 0(10%) in z* depending on 
the precise nature of the matching conditions and the 
LTB model used - even though the overall errors appear 
to be tiny. Indeed, it is not too surprising to note that 
when defining the LTB model at the centre, which ap- 
proximates the actual spacetime with radiation we are 
trying to model, there is further ambiguity: when we re- 
move the radiation, do we keep the matter density fixed 
or the curvature fixed (or the age, or h)l It sounds irrel- 
evant but this also leads to differences in (tg TB — to) of 
0(£*), and consequently tens of percent in z*. 



Lemaitre solution for one fluid 



An inhomogeneous, spherically symmetric Universe 
can be described by the Lemaitre metric [53] : 



ds 2 = _ e 2<t>(t,r) dt 2 + Q ll (t ' r) dr 2 +al(t,r)r 2 <m 2 
1 — k{t,r)r z 

(A10) 
where <Z|| = (aj_r)', ' = ^ and we will use ' = §i ■ The 
time coordinate t is the proper time for an observer at tq 
having defined cj>(t,ro) = 1. 

For a universe filled by one perfect fluid with equation 
of state p — wp, the Einstein equations can be recast 
into (see also [85] for analogous equations with different 
notation): 



Hi 



Is there any way around Eq. (A7|? For nearly homo- 
geneous radiation today, as found in [1^ (the low ft/ 111 - 1 
condition), we need the left hand side to vanish, or at 
least be much smaller than £*. Let's say we adjust our 
parameters in the LTB model such that ip™ — ^o = 0; 
then, can we adjust the matching time or any aspect 
of the void to set t^ B — t m = 0? It does not ap- 
pear so because a dust and radiation FLRW do not stay 
synchronised. Similarly, one can fine tune a matching 

redshift such that (tjj™ - t ) = (4™ - d° ut) ), (effec- 
tively, by adding an error associated with the curvature 
which cancels the error on radiation; e.g., a model with 
f2 ? n n) = 0.5, nin ut) = 0.7, h^ = 0.6 does this if wc 
choose z m « 2.7, and are careful to match the LTB scale 
factor to the FLRW one at z m \ alternatively, a model 
with flin } = 0.4, fl { 7 ° ut) = 0.7, M in ) = 0.6 matched at 
z = z m « 6.95 can achieve this), but this actually then 
implies the wrong temperatures in and out today; so this 
is no measure of accuracy. To avoid complications associ- 
ated with curvature one has to choose a matching redshift 
at z m ^> 2 — 3, as also shown in Fig. [12] In this regime, 
the error introduced by the approximated treatment of 
radiation in the matching procedure grows with redshift, 
and gives changes much larger than i* in general. 



So even though the radiation is a very small contribu- 
tion to the overall energy density at z m , it makes a vital 
difference to a model with matter and radiation when we 
integrate to high redshift. This is because we have to be 
able to synchronise the time coordinates asymptotically 
at very early times with those used at the centre at very 
late times to each other in a robust way. It appears that 
matching a dust LTB model to a radiation filled FLRW 
model is subtle and cannot be sensitive enough for the 
results we are interested in. 
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In case of pure dust, w m = and the LTB solution is 
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recovered. Indeed, Eq. (A14| implies p oc a , , so that 



the time-dependence of p is entirely encoded in a± , and 
the mass M = 8irG/3 pa\ inside a comoving shell of 
radius r is conserved. Moreover, zero (or homogeneous) 
pressure leads to <f>' — in Eq. (A16) and, by a simple 
change of time coordinate, to <f> = 0. This in turn means 
k = k(r), see Eq. ( A15 1, and the Friedmann equation can 



be rewritten in the standard way along each wordline but 
in terms of p instead of p. 

The picture is different for a radiation fluid, where 
w r = 1/3 and the evolution of p deviates from the aj 
scaling. Indeed, the inhomogeneous pressure leads to ex- 
change of entropy between neighbourhood regions and 
so to non-conservation of entropy inside comoving shells. 
Because of this, the curvature along a given worldline 
no longer stays constant as it does in an LTB universe, 
but rather its time-evolution is related to the radiation 
inhomogeneity via Eqs. (A15) and (A16). This feature 



is clearly present also in the general case of an inho- 
mogeneous universe including both matter and radiation 
components and could lead, in the matter era, to a non- 
negligible deviation from LTB dynamics depending on 
radiation profile, as we will illustrate. 
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Einstein equations including matter and 
radiation fluids 



Now we consider a universe filled by matter and radi- 
ation which have energy momentum tensors: 



rpab 
m 

i.b 



pn n 



T"" = pm a u b +^pg ab 



(A17) 
(A18) 



where p and n a are energy density and four-velocity of 
matter, p and u a are energy density and four-velocity 
of radiation, and g a b is the metric tensor. We have ne- 
glected the anisotropic pressure term for the radiation 
which couples to the hierarchy of multipole equations 
from the Liouville equation. Note that, in an inhomo- 
geneous universe, matter and radiation are not comov- 
ing, and so n a ^ u a . For the large and shallow void 
models considered here, the peculiar velocity v p between 
the two frames is however expected to be small until late 
times, namely, until the inhomogeneity crosses the hori- 
zon. We describe baryons and dark matter as a single 
matter fluid, which is an approximation since baryons 
are coupled to photons until last scattering. However, 
the peculiar velocity at LSS is typically extremely small 
and so these two matter fluids can be confidently treated 
as comoving. In the following we introduce the 3- velocity 
v a (where the peculiar velocity is v p — \Zg a bV a v b ) and 
present the Einstein equations to first order in the pecu- 
liar velocity i.e., neglecting terms 0(v 2 ) in the solution. 
(The full equations may be written down but they are 
seemingly intractable.) 

We have a choice to make for the coordinates. We 
consider two cases: a frame comoving with matter and a 
frame comoving with radiation. 
Matter frame 

The coordinates are fixed such that 

n a = (e - *, 0,0,0), 
v a = (0,v,0,0) 

u a = j(n a + v a ) : 7 = 1/v/l - v a v a , n a v a = 0. 

where v — v(t,r). In this case, V a r^ 6 = implies that 
</>' = 0, and so we can set <p — as a gauge degree of 
freedom. The equations of motion are then: 



Hi = 



2M 



k 

^2 



a\r 3 

2 1 

M = a\\a 2 ,r 2 pv pr 3 a 3 , H± 

3 II J. f* q^ ± -L 

4aiiaj_ 

k = pv 

6r 

p = -(#|| + 2H±)p-\ 
v = -(2H± -5#||)v 



2a\\a_\_r 
3(1 - kr 2 ) 
(1 - kr 2 )p' 
4a n p 



ppv 



{a\\a\r 2 p - 2M')(Ara\\a ± H ± v - 3(1 - kr 2 )) 
3a||a^_r 2 (l — kr 2 ) 



Radiation frame 

The coordinates are fixed such that 

u a = (e-*, 0,0,0), 
v a = (0,«,0,0) 

n a = 1 (u a + v a ) : 7 = 1/-/E - v a v a , n a v a = , 

(note the different definition for v) which leads to: 
2M k 



H 



a 3 , r 3 a 2 . 



M = - 



1 



a 2 ±r 2 e 4>(i _ kr 2 ) 

2^ r ~ 6 



pv - -pr ' a ± H± 



a^il-kr 2 ) a ± H ± (l - kr 2 )p' 
k = pwH 

a»r la\\rp 

4 ra ± H ± , 2a ± re* 

p = --(#|| + 2H ± )p p + — ppv 

3 3a\\ 6a\\ 

e^p' 



Ap 

9 4/i 

(a^a^p - 2M')(ra±H ± e-' t 'v - ay) 

P 2 2 O 

a y a l r 

Note the different meanings of all the variables in the 
two frames. In the matter frame the main effect of the 
radiation lies in the evolution of the spatial curvature, 
while in the radiation frame there are additional effects 
such as the de-synchronization of clocks between neigh- 
boring worldlincss directly induced by gradients in the 
radiation energy density (i.e., t is not the proper time). 
Note also that the void profile - while fixed in LTB - can 
now change in (co-moving) size. Hence, the scale of the 
radiation profile will tend to grow relative to the matter 
profile. 

In order to solve the field equations, we note that we 
have a first-order pde system in 5 variables (in the mat- 
ter frame this is a±, M, kr 2 , p, v; in the radiation frame, 
we can integrate for 4> trivially and use oj_, M, kr 2 , p, v). 
Given 5 evolution equations we must specify initial condi- 
tions for the 5 variables, as well as 5 boundary conditions. 
The important point for this paper is that the boundary 
is r = 0, and the obvious choice of boundary condition is 
v = 0,p' = 0,p' = 0. This ensures that M = (kr 2 ) = 
along r = 0, and that this worldline is FLRW. The other 
boundary is r — oo and is never affected by r — (except 
at t = oo); if the initial conditions are chosen to be flat as 
r — > oo then that will give a model which has the same 
asymptotics as we have used here. Although it might 
seem intuitive that information from the initial data can 
propagate into the central worldline and prevent FLRW 
evolution, the boundary condition there propagates out 
into the spacetime to ensure this does not happen. 

For our purposes, the area distance to LSS in the only 
quantity for which a full integration of such systems of 
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partial differential equations is needed. Indeed, as men- 
tioned many times in the paper, we compute the small- 
scale CMB by splitting the physics of decoupling and 
line-of-sight effects. For the first, equations reduce to the 
FLRW case since at LSS (r 3> void size) gradients are 
negligible (for Gaussian- like density profiles). The area 
distance cLa then accounts for line-of-sight effects. 

We motivate below that, for this quantity, LTB formu- 
las can provide a good approximation and, in the rest of 
the paper, we assumed (Ia{z) as computed in the LTB 
framework. 



4. Redshift relations for the central observer 

An incoming photon along a radial null geodesic satis- 
fies ds 2 = = dO, which leads to (in terms of the proper 
time t for an observer located at r where </>(£, r ) = 1): 



Eqs. (A22) and (A23) are general and do not rely on any 



dt 



a\\(t,r) 



^/l-k(t,r)r 2 



(t,r) dr 



(A19) 



We define \ = e~'*' a\\/vl — kr 2 . Considering another 
light ray emitted with an infinitesimal time interval r of 
delay: d(t+r) = —dr x(t+T, r), and taking the first order 
in the Taylor expansion \(t + r) = x(£) + \(t) r, one can 
easily derive dr — — x(£) r dr. Then, from the definition 
of redshift: T obs /r = l + z, it follows dr/r = —dz/(l + z) 
and (see also Ref. [5S]): 



dz 

1 + z 
dz 

l + z 



±(t,r)dr , 

x(t> r Y 



dt 



(A20) 
(A21) 



In the LTB case, xltb = ay (t, r)/\/l — k(r) r 2 and 
Xltb/ Xltb = H\\ . It implies that, along the same 
geodesic, l/a± scales roughly as l + z. Indeed 
dlog(l + z)/d\oga± — dz/dt ■ dt/da± ■ a±/(l + z) = 
— Hn iJJ , which is ~ — 1 for smooth shallow voids, where 
we defined H± = (d± + a' ± dr/dt)/a± = H± + (1 — 
a\\/ax)/(r-x)-H±. 

In the more general case considered here this relation 
becomes: 



d log (l + z) 
d log a± 



X 
XH± 
H\ 
#i 



(A22) 



1- 



Hn 2{l-kr 2 )H« 



Defining G= -H\\/H± and F = -(p/H^+kr 2 /[2 H\\ (1- 
kr 2 )], the scale factor along the past null cone of the 
observer at the centre is then given by: 



aj_(z) = cxp 



dz' 

(l + z')G(z')(l + F(z')) 



nullcone- 

(A23) 



approximation. Note that the same integral appears in 
the constraint Eq. (25). In the matter dominated era, 
G is mostly set by the matter profile. Deviations from 
— 1 are thus related to gradients in the matter density, 
so are present in the LTB scenario as well and can typ- 
ically induce a shift of 0(%) in the a± vs. z relation. 
The F-term, however, induces a deviation from the anal- 
ogous LTB relation. To quantify it, we now consider the 
solution presented in Sec. |A3| 

In the LTB coordinate frame (namely, the matter 
frame), we can set (f> = <j) = 0, which shows that the 
correction is associated to the temporal variation of the 
curvature term. This is induced by the inhomogeneous 
pressure as outlined in Sec. |A2| Combining the third 
and fourth equations of the system, it's interesting to 
note that dlogp/dt = -H\\ - 2H ± - kr 2 /[2(l - kr 2 )}. 
The latter term, which gives the correction in the matter 
density evolution with respect to the LTB computation, 
is the same as the correction in the redshift-relations. 
Using v p = an/y/l — kr 2 v, the _F-term reads: 



F 



2a±r 



3mVl-kr 2 



ix v p 



(A24) 



which is proportional to the peculiar velocity. This shows 
once again that it has its origin in the inhomogeneity 
of the radiation density and the associated pressure- 
gradient which induces a peculiar velocity in the radi- 
ation component with respect to the matter frame. 

To compute the redshift of CMB photons, we focus on 
the radiation frame, where v p = yl — kr 2 /a\\ v, and 



F 



a± r 
6 an -ffn 



H 






2e« 



pv 



3H(r )-H\\ -2H ± 
3 #|| 

(A25) 
where e^ = (/i(ro)/^) 1 ' 4 . We derive an estimate of the 
size of F by computing the evolution of p, /z, and a±,\\ 
neglecting velocity terms in the system of Sec. |A3| and 
then plugging in the solutions to derive the evolution of 
the coordinate velocity v (fifth equation) . We stress that 
such estimate only gives a rough order of magnitude idea 
about the shift in the a± vs. z relation given by inhomo- 
geneous radiation. To properly compute it, it is crucial 
to correctly estimate the evolutions of size and gradient 
of the radiation profile, which requires the integration of 
the full system of equations (strictly speaking, one cannot 
separate /i and v p evolutions). 

We focus on an example and set parameters as in the 
benchmark case described in Fig. M (with the shape of 
radiation profile the same as the matter). Results are 



shown in Fig. 14 



The main conclusion that can be drawn is that includ- 
ing inhomogeneous radiation (1 -\-z) ■ a±_ can significantly 
differ from 1 (by ~ 30% in this approximation) at z> 1. 
The correction increases corresponding to the profile gra- 
dients and then flattens in the asymptotic FLRW region. 
Critically, we need a shift of just about this amount to 
satisfy Eq. ([25 1 . 
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From another point of view this implies that we have 
freedom to choose a* where this choice means that 
selected a radiation profile giving the proper shift 
Note that it is analogous to choos- 
From the approximated picture in 



a 

we 

in Eq 

ing t* in Eq 



(A23) 
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Fig. [14] it looks plausible that a radiation profile with an 
0(1) density contrast today (similar to matter) can give 
|1 — a* • (1 + z*)\ ~ 10— 20%, which is the shift typically 
required by the fit to CMB data in the models considered 



in the rest of the paper (see Sec. Ill I and correspond to 
F ~ 1%. 

Assuming the central wordline to be FLRW, this is pro- 
vided by an inhomogeneity 0(1) in the radiation density 
at times around i* [58]. More generally, the study of 
the class of radiation models selected by a* requires a 
numerical solution of the system in Sec. A 3 (plus initial 
conditions possibly provided by a production mechanism 
for the inhomogeneity linking matter and radiation initial 
profiles) and deserves a dedicated work. We will consider 



(out) 



n (out) 



a*~ ' as a free parameter (set by T , as explained 
above). 

Similarly, a correction is expected also in the compu- 
tation of the coordinate distance r from Eq. |A20| The 
difference |1 — r/rLTs| can be sizable (0(10%) for pro- 
files and approximation illustrated here) and the shift 
goes in the opposite direction with respect to the shift in 
a* discussed above. On the other hand, to precisely state 
whether or not they cancel in the area distance <1a = a*r* 



requires the integration of the system in Sec. A3 and de- 
pends on the details of the radiation profile. It is possi- 
ble that they don't, in which case some of our results in 
Sec. |IIIB] will change; in particular, it is conceivable that 
an asymptotically flat model may be possible. We will 
simply assume cLa to follow from LTB equations, which 
is our main approximation (or, alternatively, it can be 
seen as a restriction of the analysis to a particular class 
of radiation profiles) . 

Note that, in order to estimate F, we need to com- 
pute the peculiar velocity. Our results roughly com- 
pare to Refs. [T5J [3T1 EH] 7 where the velocity is computed 
in an LTB scenario considering the dominant contribu- 
tion given by the dipole, i.e., v p /c = (AT/T)di po ie, and 
to Ref. |52j . where a covariant formalism is introduced. 
However, our result is strongly dependent on the assump- 
tion for the radiation profile (and on the approximation 
considered for the solution). Therefore we stress that 
in order to have a complete estimate of the peculiar ve- 
locity, and in turn of the kinematic Sunyaev-Zeldovich 
effect in void models, one has to take into account the 
effect of inhomogeneous radiation in geodesic equations. 
It will probably lead kSZ data to select radiation profiles 
rather than rule out matter inhomogeneity, as discussed 
in Sec. HVl 



What if the central worldline is not FLRW? 
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FIG. 14: Peculiar velocity v p (red) and |1 — ax • (1 + z)\ 
(blue) estimated in the approximation described in the text. 
The deviation of (1 + z) ■ a± from 1 is computed evaluating 
Eq. ( |A23[ ) in the radiation frame. The dashed line shows the 
result neglecting the _F-term. 



While we have argued that the central worldline can 
be considered to be FLRW to a good enough approx- 
imation (rather, we used this assumption to calculate 
early-time conditions for matter and radiation; it is not 
directly used to calculate constraints from the CMB), let 
us assume for a moment that this is incorrect. The ra- 
diation temperature at the centre is determined by the 
radiation streaming into the central worldline along null- 
cones, and so one can argue that it may be determined 
by the redshift to a surface of constant time in an LTB 
model, as opposed to the evolution of radiation along the 
central worldline. (Of course, if the full spacetime solu- 
tion were known these would be the same.) As argued 
above, we need to know precisely the surface t*(T*) to 
calculate this accurately, but the LTB approximation at 
least illustrates the key idea for this central temperature 
calculation. Let us consider where the central worldline 
evolution changes our analysis. 

In the radiation frame, the temperature always be- 
haves asToc 1/drad; where a ra d is the mean length scale 
in that frame (defined covariantly through the frame ex- 
pansion rate) . In the matter frame the temperature evo- 
lution is more complicated because of the radiation flow- 
ing through the frame. At the centre, let us assume that 
it can be calculated at any time t by T = (1 + z#)T#, 
where # denotes the surface of constant time t#, in a 
pure dust LTB model. If we do this we find, approxi- 
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mately, 



T 



To 



Too 1 



for a > a 



for a < a, 



(A26) 



where a is the value of the scale factor (at the centre) 
where the transition from normal 1/a scaling changes 
(i.e., when the central observer starts to see the void), 
and n gives the scaling behaviour at late times (n — 1 
corresponds to FLRW evolution). The consequences of 
this are that: 



The baryon-photon ratio evolves: we have tjq 



(in) 



,3(1- 



^ ln) - 



• The age calculated at the centre will be different. 



n(in 



(in) 



• Tr> and T^ 1 ' will change. 

So, if we take 7 Li constraints on rj^ in ' , for example, then 
this is the early time value, giving rji , and is less to- 
day by about a factor of two or so (depending strongly 
on the void parameters), given by a 3 ^ - ™- 1 . This cither 



requires a lower M ln ' or fl," 1 if fb — const., or a lower 
fb at the centre - see Eq. (11). In fact, if we enforce 
fb — const, and use a very low value of J?q , as esti- 
mated from this LTB scaling law, we find models similar 
to [1H1 H2J [SHI EI] • As far as the age is concerned, the 
change is tiny since radiation is subdominant when the 
change in scaling takes place (changing, e.g., fl„i by a 
fraction of a percent would adjust for this). The rest of 
our analysis goes through as discussed, and, in particular, 
the constraint Eq. (25) will be automatically satisfied. 



The fact that T* and T e q change makes no differ- 
ence to our analysis; these are determined after-the-fact 
anyway, as illustrated in Fig. [2J Of course, as we have 
discussed in this Appendix, we are unable to calculate 
T along the central worldline in the real model following 
this approach because we don't accurately know where 
t* is when we integrate from the central worldline down 
a past lightcone. Finally, we note that specifying a tem- 
perature evolution law for the central worldline is akin to 
choosing an inhomogeneous radiation profile: it is clear 
from this consideration too that we are free to do just 
this. 
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